Option Pricing from Wavelet-Filtered Financial Series 
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We perform wavelet decomposition of high frequency financial time series into large and small 
time scale components. Taking the FTSE100 index as a case study, and working with the Haar 
basis, it turns out that the small scale component defined by most (~ 99.6%) of the wavelet 
coefficients can be neglected for the purpose of option premium evaluation. The relevance of the 
hugely compressed information provided by low-pass wavelet-filtering is related to the fact that 
the non-gaussian statistical structure of the original financial time series is essentially preserved for 
expiration times which are larger than just one trading day. 
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I. INTRODUCTION 



The problem of option pricing [lH3( has been a main 
topic of investigation in much of the econophysics liter- 
ature, challenged by the well-known inadequacy of the 
standard Black-Scholes model to the real world pH^. 
Options are an imperative element in modern markets, 
since they play a fundamental role, as convincingly shown 
long ago by Black and Scholes, in reducing portfolio risk. 
As an alternative to the Black-Scholes model, one of the 
authors has implemented an option pricing scheme which 
is based on the evaluation of statistical averages taken 
over samples generated from the underlying asset log- 
return time series @. This method, which we will refer 
to as "Empirical Option Pricing" (EOP), has been suc- 
cesfully validated through a careful study of FTSE100 
options. 

A deeper understanding of the statistical features of fi- 
nancial time series is in order, since this would eventually 
allow us to replace real samples by accurate synthetic fi- 
nancial series, improving the statistical ensembles used in 
EOP. As a closely related issue, our aim in this work is to 
show that financial series can be hugely compressed (we 
mean lossy compression, in the information theoretical 
sense) by wavelet-filtering, without spoiling option pre- 
mium evaluation by EOP. The low-pass wavelet-filtered 
signal contains log-return fluctuations defined on time 
scales larger than a few hours and it is likely to yield, 
due to its high compression rate, a more suitable ground 
for modeling and synthctization. 

This paper is organized as follows. Sees. II and III 
provide brief accounts, respectively, of the EOP method 
and of the low-pass wavelet-filtering procedure that has 
been applied to our analysis of the FTSE100 index. The 
wavelet-filtered financial series, which keeps only 0.04% 
of the total number of wavelet components of the original 
signal is noted to encode the essential statistical informa- 
tion needed for a consistent evaluation of FTSE100 op- 



tion premiums with expiration times larger than a single 
trading day. In Sec. IV, we summarize our findings and 
point out directions of further research. 



II. EMPIRICAL OPTION PRICING (EOP) 

We rephrase here, without paying much attention to 
rigorous considerations, the main points of EOP @. Let 
S(t) be an arbitrary financial index modeled as a contin- 
uous stochastic process. More precisely, we write down a 
Langevin evolution equation for S(t), which is a simple 
generalization of the one underlying the Black-Scholes 
model 0,||: 



dS_ 

dt 



= fi(t)S + o-(t)ri(t)S 



(2.1) 



Above, /i(t) and a(t) are the time-dependent interest rate 
and volatility of the index S. The stochasticity of the 
financial time series comes from the gaussian white noise 
process rj(t) appearing in Eq. (|2.1[) . which satisfies to 



(v(t)) = , 

( V (t)r,(t'))=S(t-t') 



(2.2) 



Observe that both n(t) and er(t) may be regarded as 
stochastic process as well, with fluctuations correlated 
on time scales which are much larger than the correla- 
tion time of S(t). 

Working within the ltd prescription, Eq. (|2.1|) can be 
readily rewritten as 



where 

S(t) = Scxp 

with 



dt'p,(t') + x(t) - x(0) 



(2.3) 



(2.4) 
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(2.5) 
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Above, S = 5(0) is just the spot price of the index. We 
are interested, now, to evaluate the premium of an euro- 
pean option which is negotiated with strike price E and 
expiration time T . Similarly to what is done in the Black- 
Scholes treatment, where fi(t) and a(t) are constant, the 
option premium V (for, say, call options) can be obtained 
from the computation of the statistical average 

V = exp[-rT]((S(T) - E)Q(S(T) - E)) , (2.6) 

where [i(t) is replaced by r, the risk-free interest rate, in 
the definition of S(T) provided by Eqs. and (|2"3|) . 

For stochastic processes {x n } defined in discrete time, 
with time step e, like real financial time series, Eq. (|2.3[) 
can be replaced by the finite difference equation 



where 



Vn 



in 



(2.7) 



(2- 



and = ±1 is an arbitrary element of a discrete gaussian 
stochastic process, defined by (r/ n ) = and (r/ n 7] m ) = 
6 nm - From (|2.7p . we get, immediately, 



1, 

-(Xn+l 

e 



e 



and, therefore, 

1 

2 



T 1 N-l 

/ ^)^-^(fe„) 2 , 

•' Tl = 



(2.9) 



(2.10) 



where time instants are given by t n = ne, with T = Ne. 

Substituting, now, (|2.4[) and (|2.5|) (with /i(t) replaced 
by r) in (|2.6[) . we get 



S(T) = Sexp 



5 exp 



JV-l 



rT + x(T)-x(0)-±J2( Sx ^ 2 



1 



JV-l 

rT + ^ \Sx n - ^(Sx n 

n=0 



.(2.11) 



It is important to note that 5x n , which appears in the 
above expressions is, from Eq. (|2.4j) , nothing more than 
an element of the detrended log- return series, i.e., 



Sx n = ]n[S(t n+ i)/S(tn)] - ep,(t„) , 



(2.12) 



where (Sx) = 0, due to Eq. (|2~3|) . 

We are now ready to summarize EOP in the four fol- 
lowing steps: 

(i) A large period (> two years) of reasonably statisti- 
cally stationary high-frequency (minute-by-minute) log- 
return series of the underlying asset is "purified" by the 
remotion of outlier events (typically, log-return fluctua- 
tions which are larger than 10 standard deviations) and 



of the mean one-week asset's interest rate (detrending). 
The resulting series is a stochastic process {Sy n }; 

(ii) Since the historical volatility a = J ((Sy n ) 2 ) is in 
general different from the volatility of the financial series 
during the option lifetime T, we introduce a correction 
factor g to define the stochastic process {Sx n — gSy n }, 
whichyields a putative volatility a* = ga for that pe- 
riod [9J. The g- factor is the only adjustable parameter 
in EOP, which accounts for the distinction between the 
past and the future behaviors of the financial index S(t). 

(iii) An ensemble £ of samples, each of length T = 
Ne (e = 1 minute) is defined from one-hour translations 
of the initial sequence {Sxq, Sx\, 5xn-i}- In other 
words, 

£= [J {dx m A,Sxi +m A, •••j&EjV-l+mA} j (2-13) 

m,A 

where m £ N and A = 60; 

(iv) Option premiums are computed from (|2 .6[) , (|2.1ip 
and (|2.12p . with statistical averages taken over the en- 
semble £, defined in (|2.13[) . We note, furthermore, that 
the optimal value for the g-factor is found through the 
least squares method, devised for the comparison be- 
tween the market and modeled option premiums. 

A good agreement has been attained between the mar- 
ket and EOP values in a detailed study of the FTSE100 
index options @ . The comparison data is reported in the 
MKT and OP columns of Table I. 

The performance of EOP would benefit greatly from 
the use of synthetic financial series which would enlarge 
the ensemble of samples £ . Thus, one may wonder, hav- 
ing modeling aims in mind, on what are the relevant sta- 
tistical facts hidden in the financial time series. The es- 
sential question we address is, accordingly, whether the fi- 
nancial series be decomposed into relevant and irrelevant 
contributions, as far as option pricing is concerned. In 
the next section, we recall some ideas on wavelet-filtering, 
which have been crucial in the investigation of this issue. 



III. WAVELET-FILTERING AND EOP 

Log-return fluctuations are constantly affected by 
avalanches of market orders which have to do with specu- 
lative trends, and are clearly time-localized events. These 
features render the financial time series suggestively ad- 
equate for wavelet analysis. 

Since there is no requirement of continuity for the log- 
return time series, we have chosen, due to easiness of 
handling, to work with Haar wavelets [l(|. In the same 
way as for any other discrete wavelet basis, the Haar 
wavelets are labelled by two integer indices 1 < j < J 
and < k < V — 1 and are given by 



where 



^jk(t) =^ 00 (2H-k) , 



1 for < t < \ 
for \ < t < 1 



(3.1) 
(3-2) 
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FIG. 1: Monolog plots of the histograms of log-return fluc- 
tuations for both the original signal (empty circles) and the 
wavelet-filtered series (filled circles) taken for time horizons 
of (a) 100 minutes, (b) 300 minutes, and (c) 600 minutes. We 
call attention to the non-gaussian profiles of theses distribu- 
tions. 



is the function known as "mother wavelet" . Observe that 
the above basis functions are defined in the domain < 
t < 1. 



The detrended log-return series {fen, fei, fejv-i} 
of length N = 2 J+1 and zero mean [Tl| can always be 
expanded in wavelet modes as 



,7 2 J -1 



(3.3) 



j=0 k=0 



Low-pass wavelet-filtering can be straightforwardly im- 
plemented from the expansion ()3.3|) by retaining the 
modes which have the scale index j < j*, where j* is 
an arbitrarily fixed threshold. We have taken (following 
the prescriptions given in step (i) of EOP, as discussed in 
Sec. II) a financial time series of 241664 minutes (around 
two years of data) for the the FTSE100 index, ending on 
17th novembcr, 2005. The series is partitioned into 59 
subseries, each of length 4096 (corresponding to J = 11 
and about two weeks of market activity) , which are then 
wavelet-filtered with threshold parameter j* — 4 (com- 
pression rate of 99.6%). 

Since wavelets with j > have zero mean, we expect 
that the histograms of log-returns ^2^=o w ' n not ^ c 
much affected for time horizons T = Ne > 4096/2 : ' = 
256 minutes. This is actually verified in Fig. 1. There- 
fore, on the grounds of Eq. (|2.11[) . it is clear that option 
prices can be alternatively estimated through the use of 
the low-pass wavelet filtered scries within EOP for time 
horizons which arc larger than one trading day (T = 510 
minutes). 

In Table I, we report the computed premiums for call 
options based on the FTSE100 index with expiration 
times ranging from a few days to one month, in deccm- 
ber 2005 and january 2006. The agreement between the 
original and the wavelet-filtered option premium evalua- 
tions is significative. It is important to recall, as already 
indicated in Ref. @, that the Black-Scholes framework 
is unable to yield good estimates of the market option 
premiums listed in Table I. 



IV. CONCLUSIONS 

We have found, taking the FTSE100 index as a case 
study, that its high frequency (minute-by-minute) time 
series can be highly compressed for the purpose of option 
pricing. The original and the low-pass wavclct-filtcrcd sc- 
ries have remarkably similar performances within EOP, 
even for a compression rate of 99.6%, which means that 
only 967 out of the original 241664 wavelet components 
have been selected throught the wavelet-filtering proce- 
dure. The retained wavelet coefficients have scale index j 
smaller than the fixed threshold j* = 4 and are associated 
to log-return fluctuations defined on time scales larger 
than a few hours. It turns out, thus, that one is entitled 
to use the filtered time series to precify FTSE100 options 
with expiration times which are larger than just one trad- 
ing day where log-returns are still clearly non-gaussian 
random variables. A promising approach to option pric- 
ing, deserved for further investigation, is to address the 
problem of series synthetization from the analysis of the 
statistical properties of the compressed wavclct-filtcrcd 
financial indices directly in wavelet space, in a spirit simi- 
lar to what is done in the context of artificial multrifractal 
series [Hj]. It is likely that EOP would, then, be consid- 
erably improved from the use of much larger synthetic 
statistical ensembles. 
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TABLE I: Call option premiums taken from the market (MKT) are listed together with the EOP evaluations performed with 
the original (OP) and the wavelet-filtered (OP) series. The market data was recorded on 02dec05, 06dec05, 09dec05, 19dec05, 
03dec06, and 12dec06; spot prices are indicated by S; the respective g-factors (and putative volatilities; see Sec. II) are 
g = 0.81 (a* = 6.1%), g = 0.91 (a* = 6.9%), g = 0.94 (a* = 7.1%), g = 0.78 (a* = 5.9%), g = 0.83 (a* = 6.3%), and g = 0.78 
{a* = 5.9%). The risk-free interest rate is r = 4.5%. The three first dates in december refer to options which expired on 
16dec05. For the other three dates, expiration date is 20jan06. The mean volatilities measured between 02dec05 and 09dec05 
and between 19dec05 and 12jan06 were a = 8.0% and a = 6.1%, respectively. The brackets in some of the OPo evaluations 
indicate option premiums which are probably affected by poor sampling of the underlying financial series. 
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